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Abstract 

An approach for describing the hindrance of the nuclear 2v(3(3 -decay amplitude is 
proposed. The approach is based on a new formula obtained by a model-independent 
transformation of the initial expression for the amplitude. This formula takes ex- 
plicitly into account the hindrance of the decay-amplitude due to the presence of the 
collective Gamow- Teller state. Calculations are performed within the simplest ver- 
sion of the approach. Calculated and experimental 2v(3f3 halve-lives are compared 
for a wide range of nuclei. 
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1. Introduction 



A rough analysis of a great body of 2vf3f3 -decay data shows that the nu- 
clear decay-amplitude is hindered as compared to the one evaluated within the 
independent-quasiparticle approximation, (see e.g. refs [1,2]). The hindrance 
is caused by the presence of the isobaric analog and Gamow- Teller nuclear 
collective states (IAS and GTS) at a rather large excitation energy. These 
states exhaust a major fraction of, respectively, the Fermi and Gamow- Teller 
strength. To describe the hindrance, different versions of the quasiparticle 
random phase approximation (QRPA) are used in calculations of the nuclear 
2vf3f3 -decay amplitude. It has been found (see e.g. refs.[3]-[6]) that the calcu- 
lated amplitudes are very sensitive to the particle-particle interaction strength. 
For this reason the predictive power of the current QRPA versions is poor. 
A similar instability of the calculated OvfiP - decay amplitude has also been 
found [7,8]. 

To bypass this difficulty, which probably is caused by an inconsistency in the 
current versions of QRPA, we propose a rather different approach to describe 
the hindrance of the nuclear 2v(3j3 - decay amplitude. The method is based 
on a new formula obtained by a model-independent transformation of the ini- 
tial expression for the amplitude. This formula takes explicitly into account 
the decay-amplitude hindrance due to the presence of the collective GTS. The 
approach is similar to the one widely used in the description of Fermi-type 
excitations in intermediate and heavy mass nuclei. These excitations can be 
described with the explicit use of the experimental fact, that the IAS exhausts 
most (more than 95%) of the Fermi sum rule N — Z. For these nuclei the mean 
nuclear Coulomb field is mainly responsible for the redistribution of the Fermi- 
strength from the IAS to its satellites. If the variable part of this field were 
neglected, the IAS would almost exhaust 100% of the Fermi strength and 
the Fermi 2v(3(3 -amplitude would vanish for transitions to the ground and 
low-excited states of the product nucleus. This approach, which allows one to 
evaluate the small Fermi-strength of the IAS satellites populated virtually in 
the 2vf3f3 -decay process, has been originally proposed by Lane [9] and is 
based on the approximate conservation of isospin symmetry in nuclei. 
A similar situation occurs in the description of Gamow- Teller excitations. It 
is experimentally established that in intermediate and heavy mass nuclei the 
GTS exhausts about 70% - 80% of the sum rule 3ejj(7V — Z), where e q ~ 0.8 
is "an effective charge" describing, in particular, the re- normalization of the 
axial coupling constant for the weak interaction in nuclei. The GTS collec- 
tivity, directly related to this experimental fact, is explicitly used in the pre- 
sented work to describe the hindrance of the nuclear 2vf3(3 -decay amplitude. 
Early attempts along this line have been undertaken in refs. [10]- [12] within 
a perturbation theory. As a first step in applying the proposed approach, the 
decay amplitudes are calculated by the new formula using the BCS and the 
pair-vibration models to describe the corresponding states of isobaric nuclei. 
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Calculations performed for a wide range of nuclei show clearly the hindrance of 
the decay amplitude due to the presence of the collective GTS. The 2vf3~ f3~ 
and 2b / f3 + j3 + (including the electron capture) half-lives are calculated using 
model parameters taken from independent data. Results of calculations are 
compared with available experimental data. 



2. The basic formula for the nuclear 2v(3(3 - decay amplitude 



To make the following derivation of the new formula for the 2v(3(3 -decay 
amplitude more transparent, we will use the same approach for the Fermi and 
the Gamow- Teller transitions. Let G^' = J2 a 9^ be the operators of allowed 
Fermi (g^ = r^) and GT (g^ = ar^) (3 - transitions. First, we consider 
the 2u(3~(3~ -decay amplitude. Let \i) and |/) be the ground state wave 
function, respectively, of a double-even parent nucleus (Z, N + 2) and of a 
product nucleus [Z + 2,7V), and Ef are the corresponding energies. The 
expression for the nuclear 2v(3(3 - decay amplitude has the form (see e.g. refs. 
[13,14]): 

M G = Y J {f\G { - ) \S){S\G^\i)u s \ (1) 
s 



where the sum is taken over intermediate states of an isobaric nucleus (N + 
1, Z + l) and ujs = E s — (Ei + Ef)/2 is the excitation energy of an intermediate 
state. The 2u(3 + (3 + - decay amplitude is equal to conjugate value Mq. 

We transform eq.(l) in the following way. Let \G) be one of the intermediate 
states with the energy Eg and V G ^ be the following operator: 

V^"> = -ifcG(-> - A G G*-\ A G = E G -E i , (2) 



where a dot over the operator implies the time derivative. Matrix elements of 
this operator can readily be constructed: 



-(u G -ujs)(S\G^\i) = (S\V^ ) \i) , 

-(uj G + uJs)(f\G^\S} = (f\V G - ) \S) , (3) 

giving (G\V G = 0. Using these matrix elements we transform amplitude 
(1) to the expression: 

m g= 1^ — —77-2 — ~2\ ^2 • yv 
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The last term in this expression can be transformed with the help of eqs.(2),(3) 
and the equalities 

= (f\[vt\G^]\i) = E((/I^V)(^|G ( - ) K) - (f\G^\S)(S\vt ] ^).(5) 

s 

The commutator vanishes because of the fact that (r^) 2 = 0. As a result 
of the transformations performed with the help of eqs.(2)-(5) we get the final 
expression for the nuclear 2uf3f3 - decay amplitude: 

M G = ^Y,(f\yt ) \S)(S\vt ) ^ s 1 . (6) 

s 

Summation in this equation formally includes the state \G) as well, because 
(G\vt ) \t} = 0. 

Eqs. (1) and (6) are equivalent. However, the 2vf3f3 -decay hindrance due to 
the presence of the collective states (IAS and GTS) can be explicitly taken 
into account in eq.(6) provided that the state \G) is, respectively, taken as 
one of these states. Indeed, if the above states each exhaust 100% of the 
corresponding particle-hole strength, i.e. \G) ~ G^\i), then matrix elements 
(S\V G ~^\i) and, therefore, the amplitudes M G are equal to zero. It is now clear, 
in a model-independent way, that if the collective state \G) exhausts most of 
the sum-rule strength, the 2vf3f3 -decay will be strongly suppressed. 

For the above reasons the use of eq.(6) allows one to bypass the difficulty con- 
nected with the use of current versions of QRPA which are based on eq.(l). To 
show it we first consider nuclei without nucleon pairing using the phenomeno- 
logical nuclear mean field U and the isovector part of the Landau-Migdal 
particle-hole interaction F. The mean field contains isoscalar, isovector, spin- 
orbit and Coulomb parts: 

U = U (r) + U T (r)r^ + U al (r)(al) + U c (r)(l - t®)/2 . (7) 
The interaction 

F=(F T + F^a^mSin - r 2 ) (8) 

contains two phenomenological parameters F g (g = t,<tt) which are widely 
used within in the theory of finite Fermi-systems [15]. Let p(r) = (p^(r) + 
r( 3 )p( - )(r))/2 be nucleon density: p^(r) = p n {r) ± /^(r), p n and p p are the 
neutron and proton densities, respectively. The following conditions 

U T (r) = F T p(-); U c (r) = e 2 J ff( ri )\r - r^dr, (9) 
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allow one to calculate two terms in eq.(7) in a self-consistent way (the first 
condition is discussed in e.g. ref.[16]). 

Let us turn to the operator = £ a Vj-\a) (2). It can be shown using the 
coordinate representation of the RPA equations that: 

-ihg^ = [h,g^]-2F g [p,g^\ , 
Vj"> = [h,g^\ - 2F g [p,g^\ - A G g^\ (10) 

Here, h = t + U is the single-particle Hamiltonian. Assuming \G) ~ G^\i), 
we find an approximate expression for the energy A G using the equality 
(G\V G ~^\i) = and eq.(2): 

a ~ m l+ \G^]\i) 

Ag = - lh (i\[GW,G(-)]\i) ■ 

This expression can also be obtained by means of the energy weighted sum rule 
where Ac equals to the mean energy of Fermi or GT excitations, respectively. 
This formula has acceptable accuracy because the IAS and GTS exhaust most 
of the corresponding particle-hole strength. 

For the case of Fermi transitions (g^ = t^) we obtain from eqs.(7)-(ll): 

= (U c (r) - A c )r(-); A c = (N - Z + 2)~ l J U c {r)p^\r) df .(12) 

These well-known equations are a result of the approximate conservation of 
isospin symmetry in medium and heavy nuclei, where the mean Coulomb 
field is the main source of violation of this symmetry. Similar calculations 
for the case of Gamow- Teller transitions (g^ = gt^) lead to the following 
expression: 

= (U c (r) - A o )0V<-> + (U al (r)[al,a\ - A al a)r^ + 

{25F P { -\r)-A 5F )aT { -\ (13) 

where A c is determined by eqs.(12), 5F = F T — F ar , 



A al = ~(N-Z + 2)- 1 (*| £ U al (r a )aJ a \i), 

a 

A SF = 2{N-Z + 2)~ 1 5F J (p ( "Hr)) 2 df (14) 
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and 



A aT = A C + Ki + & SF 



(15) 



The three terms in eq.(13) arise as a result of spin-isospin symmetry violation 
(this symmetry may be called a "projection on the charge-exchange channel" 
of the spin-isospin SU(4) - symmetry [17]). However, due to the smooth radial 
dependence of both the mean Coulomb field and the excess neutron density the 
spin-orbit part of the nuclear mean field is the main source of the redistribution 
of the Gamow- Teller strength from the GTS to its satellites. In other words, 
only the second term in eq.(2) needs to be taken into consideration in the 
analysis of the Gamow- Teller 2vf3f3 - decay amplitude by means of eq.(6). 
The other terms are only important for the calculation of the energy of the 
GTS (15). 

The equations equivalent to eqs. (10) have not yet been formulated for nuclei 
with pairing. Nevertheless, one can state that the operator (12) is not 
changed since paring forces conserve isospin symmetry. In the case of Gamow- 
Teller transitions the operator (13) is expected to be somewhat modified 
due to that part of the particle-particle interaction, which violates the "pro- 
jection" of SU(4) - symmetry. The relative contribution of this interaction to 
Vj^ can be roughly estimated as the ratio of the energy gap to the mean 
single-particle spin-orbit energy-splitting, which is rather small (no more than 
20%). This contribution is omitted in the following analysis. 

In view of a higher degree of the isospin-symmetry conservation as compared 
with the spin-isospin-symmetry conservation one can conclude that Mqt ^ 
Mp. For this reason we consider in the following only the Gamow- Teller 2z//3/3 
- decay amplitude. 



3. Evaluation of the decay amplitude within a simplest version of 
the approach 

As a first step in evaluating the amplitude M GT with the help of the eq. (6) 
we use the simplest approximation for describing the corresponding states in 
isobaric nuclei. Namely, we use the BCS model for describing subsystems with 
strong nucleon pairing and the pair- vibration model for describing "magic ± 
two nucleons" subsystems. 

Within this approximation the nucleon densities (r) are equal to 

P {± \r) = i- (lX(0(2> + IK ± IX(0(2*r + IK) , (16) 
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where R\(r) are the radial single-neutron (A = v) and single-proton (A = ir) 
wave functions; A = n r , l,j is the set of the single-particle quantum numbers, 
and n x are occupation numbers satisfying the equations: 

X)(2j„ + lK = J /V + 2, ^(2j w + lK = Z. (17) 

For nuclei with strong pairing in any nucleon subsystem the occupation num- 
bers are n x — v x — 1 — u x , where v x and u x are the Bogoliubov-transformation 
coefficients. In the case that a nucleon subsystem is the "magic ± two nucle- 
ons" subsystem the occupation factors are: n x — n™ + c\{l — 2n™), where n™ 
are the occupations numbers for the magic subsystem; c x are the coefficients, 
which determine the pair-vibration wave function and satisfy the normaliza- 
tion condition: £a c|(1 - 2n 1 x n ){2j x + 1) = ±2. 

The terms A c , A CT j and A^ determining the GTS energy (15) can also be 
calculated easily according to eqs.(12),(14). In particular, we have 

A CTi = ~(N -Z + 2)- 1 n x (\\U al (r)\\)( 3x - l x )l x (l x + 1), (18) 

X=n,u 

where (ir\U(ri(r)\v) = J R n (r) R u (r)U a i(r)r 2 dr and n x are occupation numbers 
satisfying eqs.(17). 

Calculation of the amplitude Mqt according to eq.(6) within the framework 
of the BCS model results in the expression: 



M GT = e^Ts£ (2^ + l)ti*-J»)(*\U S o(r)\v) - A so (n\u) 

TT,U \ J 

x(7r||o-||z/) 2 M 7r t; 7r M^w~J . (19) 

Here, (7r||er||i/) is the reduced matrix element; u nv = £ w + E v is the excitation 
energy of the two-quasiparticle state, £ x = \J (e x — jj) 2 + A 2 is the single- 
quasiparticle energy for the subsystem with strong nucleon pairing, e x is the 
energy of single-particle level, and A is the energy gap. If the nucleon sub- 
system is magic in final (initial) state and "magic ± two nucleon" in initial 
(final) state, we use c x instead of u x v x in eq. (19) and S x = |A + e x — e±\ for 
particle pair- vibrations or £ x = |A + e — e x \ for hole pair- vibrations, where 
ei (eo) is the energy of the first empty (last filled) single-particle level in the 
magic subsystem. Here, 2A is the pair-vibration state energy, which coincides 
with the pairing energy P for the subsystem "magic ± two nucleons". The 
hindrance factor her can be estimated as follows: her = M gt /Mq T , where 
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amplitude 

M GT = e l ^2{^W) 2 {^\W\W) 2 UnVnUuV,/V~^ . (20) 

is evaluated according to the initial equation (1) with the use of the BCS (or 
BCS and pair- vibration) model without consideration of the hindrance caused 
by the presense of the collective state. 



4. Results of the calculations and summary 

The parameterization of the isoscalar part of the phenomenological nuclear 
mean field including the spin-orbit term is given e.g. in ref. [18]. The strength 
F T = 300 MeV fm 3 in eq. (8) is chosen to describe the experimental neu- 
tron and proton binding-energy difference for nuclei with rather large neu- 
tron excess 48 Ca, 68 Ni, 132 Sn, 208 Pb [19], for which the difference is mainly 
determined by the symmetry potential and mean Coulomb field. The strength 
F aT = 255 MeV fm 3 in eq.(8) is chosen following the suggestion of ref. [15]. 
The parameters /x n;P and A np , obeying the constraints of eqs. (17), were found 
for each subsystem with strong nucleon pairing to describe the experimen- 
tal pairing energies P = 2minA^A taken according to ref. [19]. The pair- 
vibration state energies 2A = P are calculated using experimental pairing 
energies [19]. Coefficients c\ determining the pair- vibration state wave func- 
tion for a nucleon subsystem "magic + or - two nucleons", are calculated 
using c x ~ (1 - 2n x n )/(e x - e x + 2 A) or c A ~ (1 - 2n^)/(e - e A + 2A) with 
taking into account the normalization conditions given above [20] . For all con- 
sidered nuclei the single-particle basis including all bound states as well as the 
quasi-bound states up to ~ 5MeV is used. If the nucleon subsystem in a par- 
ent nucleus is "magic ± two nucleons", and in the product-nucleus is "magic 
± four nucleons", the calculations were performed within the framework of 
the BCS model with the use of u\,v\ factors calculated for a "magic ± four 
nucleons" subsystem. 

Results of the calculations and corresponding experimental data are given in 
Table 1. The Mqt amplitudes calculated using eq.(19) are given in column 

5. The Mq T amplitudes were calculated from eq. (20) with the use of the 
same model parameters and the same number of basis states. The calculated 
hindrance factors Yiqt are given in column 6. Half-lives Tyif calculated by 
formula (Ti/2)^ 1 = G2v\Mgt\ 2 are given in column 8. The lepton factors 
taken from refs. [2,13,14] are also given (column 3). We calculated also the 
2vf3 + [3 + half-lives (including the electron capture) for those nuclei, for which 
relevant experimental data is expected to be forthcoming. 
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In this work the initial formula for the nuclear 2v(3(3 - decay amplitude 
is transformed in a model-independent way to an expression, which takes 
explicitly into account the decay-rate hindrance due to the presence of the 
collective GTS. To apply the transformed formula in its simplest version we 
use: (i) the fact that the GTS exhausts the most of the Gamow- Teller strength; 
(ii) the BCS and the pair-vibration models for describing the corresponding 
states of isobaric nuclei. The formula for the amplitude Mqt, obtained with 
the use of the above approximations , is the sum of the positive-sign terms and, 
therefore, is stable to reasonable variations of model parameters. Calculations 
are performed for a wide range of nuclei with the use of the parameters found 
from independent data. Except for Te isotopes the calculated amplitudes Mqt 
within the factor 2-3 are in agreement with the corresponding amplitudes 
Mqt deduced from experimental data. 

The present approach for evaluating the nuclear 2vf3f3 -decay amplitude 
is incorporating in a model-independent way the mechanism for the rather 
strong suppression of the decay-rate and is thus to be preferred over conven- 
tional approaches. The use of the present approach in QRPA requires further 
investigation. 
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Table 1. Calculated Mqt and Ti/ 2 values in comparison with the correspond- 
ing experimental data. 150 Nd is considered as a spherical nucleus. The lepton 
factors G 2v and calculated hindrance factors are also given. m e is the electron 
mass. 



parent 
nucleus 


type of 
decay 


G 2v 
years -1 m 2 




M GT 
m~ x 


her 


rpexp. 
1 l/2 

years 


n^calc. 
J l/2 

years 


76 Ge 


(3-(3- 


1.317 


X 


io- 


-19 


[13] 


0.0919 


0.0387 


0.131 


0.9 x 10 21 


[21] 


5.0 x 


10 21 
















0.0737 






1.43 x 10 21 


[22] 






78 Rr 


ec ec 


1.957 


X 


10" 


-21 


[14] 




0.0285 


0.090 






6.2 x 


10 23 




(3 + ec 


1.174 


X 


10" 


-21 


[14] 












1.0 x 


10 24 


82 Se 


P~P~ 


4.393 


X 


10" 


-18 


[13] 


0.0459 


0.0295 


0.093 


1.08 x 10 20 


[23] 


2.6 x 


1Q 20 


96 Zr 


P~P~ 


1.953 


X 


10" 


-17 




0.0362 


0.0678 


0.324 


3.9 x 10 19 




1.1 x 


10 19 


96 Ru 


ec ec 


6.936 


X 


10" 


-21 


[14] 




0.1005 


0.338 






1.4 x 


10 22 




(3 + ec 


1.148 


X 


10" 


-21 


[14] 








> 6.7 x 10 16 


[24] 


8.6 x 


10 22 


100 Mo 


P-P- 


9.553 


X 


10" 


-18 


[13] 


0.0954 


0.1606 


0.329 


1.15 x 10 19 


[25] 


4.1 x 


10 18 


l06 Cd 


ec ec 


1.573 


X 


10" 


-20 


[14] 




0.1947 


0.319 






1.7 x 


10 21 




f3 + ec 


1.970 


X 


10" 


-21 


[14] 








> 6.6 x 10 18 


[26] 


1.3 x 


10 22 


ll6 Cd 


P-P- 


8.000 


X 


10" 


-18 


[2] 


0.0754 


0.0788 


0.258 


2.25 x 10 19 


[2] 


1.2 x 


10 19 


124 Xe 


ec ec 


5.101 


X 


10" 


-20 


[14] 




0.0528 


0.085 






7.0 x 


10 21 




(3 + ec 


4.353 


X 


10" 


-21 


[14] 












8.2 x 


10 22 


128 Te 


p-p- 


8.624 


X 


10" 


-22 


[13] 


0.0123 


0.0529 


0.090 


7.7 x 10 24 


[27] 


4.1 x 


10 23 


i30 Te 


P-P- 


4.849 


X 


10" 


-18 


[13] 


0.0087 


0.0468 


0.085 


2.7 x 10 21 


[27] 


9.4 x 


10 19 


130 Ba 


ec ec 


4.134 


X 


10" 


-20 


[14] 




0.0568 


0.082 


> 4 x 10 21 


[28] 


7.5 x 


10 21 




(3 + ec 


1.387 


X 


io- 


-21 


[14] 








> 4 x 10 21 


[28] 


2.2 x 


10 23 


136 Xe 


P-P- 


4.870 


X 


10" 


-18 


[13] 


0.0299 


0.0341 


0.088 


> 2.3 x 10 20 


[29] 


1.1 x 


1Q 20 


136 Ce 


ec ec 


3.988 


X 


10" 


-20 


[14] 




0.0512 


0.081 






9.6 x 


10 21 




f3 + ec 


6.399 


X 


10" 


-22 


[14] 












6.0 x 


10 23 


150 Nd 


P-P- 


1.200 


X 


10" 


-16 


[13] 


0.0221 


0.0642 


0.182 


1.7 x 10 19 


[30] 


2.0 x 


10 18 
















0.0304 






9 x 10 18 


[31] 







The table should be put in the end of Sect. 4 
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